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We study hard dimers on dynamical lattices in arbitrary dimensions using a random tensor model. 
The set of lattices corresponds to triangulations of the d-sphere and is selected by the large N limit. 
For small enough dimer activities, the critical behavior of the continuum limit is the one of pure 
random lattices. We find a negative critical activity where the universality class is changed as 
dimers become critical, in a very similar way hard dimers exhibit a Yang-Lee singularity on planar 
dynamical graphs. Critical exponents are calculated exactly. An alternative description as a system 
of "color-sensitive hard-core dimers" on random branched polymers is provided. 
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I. INTRODUCTION 

Statistical mechanics on random, or dynamical, lattices consists in averaging a lattice model over some ensemble 
of graphs. One usually works in a grand-canonical ensemble, where the lattice chemical potential is the weight for 
adding a lattice site. The thermodynamic limit is the regime where the partition function is dominated by graphs 
with a huge number of sites, i.e. when the expectation value of the number of sites goes to infinity. In this regime, 
the partition function develops some singular behavior and one says that the lattice becomes critical [l], Q • 

Hard dimers are objects which can be attached to links of a lattice with some exclusion ride. Together with Ising 
spins 0] jthey are one of the first systems of statistical mechanics which have been studied on random two-dimensional 
lattices [4|. Using a matrix model formulation, the partition function can be calculated exactly. In the continuum 
limit, dimers and lattices can become critical together, and this gives rise to the Yang-Lee singularity, a non-unitary 
conformal field theory (CFT) with central charge c = —22/5 coupled to two-dimensional quantum gravity This 
universality class is generically observed in models of hard objects 0. Thus, hard dimers provide a typical model for 
multicritical behaviors on random 2d surfaces, which has paved the way to higher order multicritical behaviors H -fll 



Random rank-d tensor models have been introduced to generate random lattices in dimensions d > 2 [12Hl4j . A 
class of such models which relies on the use of an ensemble of colored graphs |15Ml7| , is now known to have a continuum 
limit in the large N scaling limit, N being the size of each tensor entry, [1 81 ] - Beyond the criticality of pure random 
lattices, other universality classes have been found [3 HO], but little has been said about their physical origin. We 
also refer to [2l[ for the appearance of non-universal behaviors, [22j for a detailed review of tensor models and [2(| 
for a synthetic presentation of the status of the field based on the result of [23| . 

To understand the physics of these multicritical behaviors, it was proved fruitful in two dimensions to consider the 
Ising model and hard dimers. Using a two-matrix model, it was found that Ising spins in the continuum have the 
universal behavior of the unitary Ising CFT (central charge c = 1/2) coupled to quantum gravity d, [24|. However in 
higher dimensions, using the large TV limit of a colored tensor model, we found that there is no phase transition at 
finite temperature [25| in agreement with previous numerical simulations (26|. 

The study of dimers coupled to random tensors has been initiated in 27j. There the model features dimers 
with a modified exclusion principle (they are not completely hard-core) and exhibits in a remarkably simple fashion 
multicritical behaviors. In particular an example of a phase transition was proposed at d = 6. Dimers have also 
been recently considered to reach multicritical points in two-dimensional causal dynamical triangulations (28j (see 
also H|). 

In the present paper we continue this investigation with truly hard-core dimers. In Sec. [II] we review the two- 
dimensional case and present the generalization to arbitrary dimensions in terms of tensors. The model is solved to 
leading order in the large N limit in Sec. IIIIl While the analysis holds for any d we have chosen to give systematic 
details on the case d = 3. Our result is similar to the two-dimensional situation: dimers become critical for some 
negative activity, sending the continuum limit to a universality class different from pure random lattices. This 
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universality class is the one already observed in [20j, |27|. Its entropy exponent 7, which characterizes the proliferation 
of microstates (graphs dressed with dimers) in the continuum limit, is the case to = 3 of the series 7 = 1 — 1/m of 

mm. 

The dominant graphs of our model are known to be in one-to-one correspondence with branched polymers (l8j . This 
provides an alternative representation of the model where hard-core dimers are mapped to color-sensitive hard-core 
dimers, as we explain in Sec. IIVI 

In Sec. [V] we revisit the model of 27j through a different generating function which is more natural from the 
view of the thermodynamic limit. The qualitative conclusions are the same and so are the exact entropy exponents 
7 = 1 — 1/m, which are also the same as for multicritical branched polymers [30j . We discuss the case d = 3 in details. 



II. HARD DIMERS ON RANDOM SPHERICAL LATTICES 



A dimer on a lattice G is a bond on a link between two nodes. It has an activity z which is the weight for putting 
a dimer on a link. In addition, one often considers some exclusion rules, i.e. constraints preventing the dimcrization 
of all lattice links. These exclusion rules create some frustration and there are usually many equivalent ways to put a 
fixed number of dimers on the lattice. A dimer configuration is a set of lattice links carrying dimers, not necessarily a 
covering. Denote T>(G) the set of configurations on G allowed by the exclusion rules. Then the partition function is 

E G (z)= J2 * W > (2- 1 ) 

DeV(G) 

where \D\ is the number of dimers in the configuration D. Making the lattice dynamical means that we now sum over 
a set of graphs, with a fixed number of vertices, 

S„(z)= Sg„(*)- (2.2) 

{G„} 

Introducing some chemical potential g for the number of lattice sites we get a grand-canonical ensemble 

00 

E(z,g) = J2^n(z)9 n . (2.3) 

n=0 

Assuming this generating function has a finite radius of convergence g c {z), depending on the activity z, the thermo- 
dynamic limit is the regime where analyticity is lost, so that the grand-canonical partition function is dominated by 
graphs with an infinite number of sites. Close to g c {z), the expectation value for the number of sites calculated in the 
grand canonical ensemble goes like (n) ~ l/\g — g c (z)\. 

The large n asymptotics of S„(z) is expected to be of the form S ra (z) ~ An J ~ 3 [g c (z)]~ n where 7 is the entropy 
exponent. In the grand-canonical ensemble, it translates to the following singular behavior S(z, g) ~ (g — g c (z)) 2 ~ 7 . 



A. The 2d case: Hard dimers on random planar graphs 

Such generating functions on two-dimensional graphs can be computed as the free energy of some hermitian matrix 
models. Taking the graphs in {G n } to be 4-valent graphs on the 2-sphere, hence planar, we have Q 



,N 2 S( Z ,g) 



dAdB exp 



-JV tr -A 2 - -A A + -B 2 - g^/z B A 3 
2 4 2 



(2.4) 



at the leading order for large N. Using the saddle point equation for the spectral density, one finds that for most 
values of the dimer activity z the continuum limit lies in the universality class of pure random lattices (i.e. pure 2d 
quantum gravity). The singular behavior of S(z,p) reads 



2 s ing(z,5) ~ (g-g c (z)) 2 7 , with 7 = -1/2. 



(2.5) 



However, when the activity reaches a critical value z c = —1/10, dimers become critical too and change the universality 
class, 



-single, g) ~ (g - g c (z c )) 2 7 , with 7 = -1/3. 



(2.6) 
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When one approaches the singularity along the curve of the continuum limit g c (z), another critical exponent is found 

^-\ogg c {z)~{z-z c ) x l 2 . (2.7) 
dz 

This critical point is identified with the Lee- Yang CFT, of central charge c = —22/5, coupled to Liouville gravity 
j HI El- A universality check can be done, using planar triangulations instead of quadrangulations, leading to the 
same critical exponents. 



B. Dynamical spherical triangulations in arbitrary dimensions via random tensors 

We now propose to generalize this approach to lattices in higher dimensions, using random tensors instead of 
random matrices. The most naive generalization of matrices to tensors, introduced in 1991, is not known to admit a 
large N expansion. Instead we need a slightly refined framework which docs support a large ./V expansion dominated 
by some triangulations of the d-sphere, known as colored tensor models. 

We introduce a collection of (d + 1) tensors (Ti*... a , i )i=o 1 ... 1 rf 5 where (i) is called the color index and the indices 
ai,...,a<j are the usual tensor indices, ranging from 1 to N. These tensors are complex and their conjugate are 
denoted f( % K These tensors are useful to generate bipartite (d + l)-colored graphs. Lines of a graph are generated by 
propagators between some T« and and hence carry a color index. This leads to the quadratic part of the tensor 
action 

d 

Squad = y ] y ] T ai ... ad T ai ... ad . (2-8) 
i— ai,...,ad 

The sites of each graph arc chosen of degree (d+ 1) with distinct adjacent colors. A set of 'black' vertices is generated 
by 

d 

S+ = Vff ^2 Hotiii— n 4i+ D (2-9) 

and a set of 'white' vertices by 

d 

{TK f }i=0 

Colored graphs obviously possess vertices and lines but also a richer structure due to the coloring. One can for example 
define faces as the maximally connected subgraphs with exactly two colors. These faces are interesting because, just 
like in random matrix models, the amplitude of a graph has a free sum per face, hence part of its scaling going like 
N F , where F is the number of faces. 

To have a well-defined large N limit, one must compensate in some way the above scaling with the number of faces. 
To this aim, it is necessary to re-scale the action with some power of N, which helps peak the integral appropriately, 
as follows 

e -N*F(g) = J ^ T (i) df (^ exp _ N d/2 + S++S _^j (2.H) 

Then the amplitude of a connected graph scales like N F 1 ~tt^ v ; where V is the number of vertices. It can be shown 
that F — IMzll y i s bounded by d and writes F — d \ d ~ x > V = d — jj^j^uj(G), where U){G) is a positive integer called 
the degree of the graph G. This provides the free energy of the model F(g) with a large N expansion 

F( g )= n-tA*« U9) = Mg) + ^f u= o^i(g) + --- ■ (2.12) 

Each / w (<?) is the free energy at fixed degree ui. For our purpose it is not necessary to give the definition of the degree 
of a colored graph, as we will only focus on graphs of vanishing degree and provide a precise characterization of them. 
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Figure 1. A typical melonic graph for d — 3 with two external legs. Notice that vertices come by pair (a black and a white, 
connected by d lines with melonic insertions). 



Note however that this framework makes sense for any d > 2. When d = 2, i.e. tensors reduce to matrices, the degree 
ui is exactly the genus of the graph, as expected. 

As in 2d, the colored graphs generated by such colored tensor models have a natural topological interpretation in 
dimension d. To each vertex of degree (d + 1), one associates a simplex of dimension d. It has (d + 1) simplices 
on its boundary, which are represented on the graph by the lines meeting at the vertex. Hence, these boundary 
simplices are identified by a color (i). A (d — 2)-dimensional sub-simplex of a <i-simplex is shared by exactly two 
(d — l)-simplices on the boundary, say of colors (i) and (j), and is therefore uniquely identified by the pair of colors 
(ij). This identification scheme of sub-simpliccs in a simplex works for any dimension and is made possible thanks 
to colors. When a line of the graph joins two vertices, this means that two d-simplices are glued by identifying one 
of their boundary (d — l)-simplex. This simplex is again uniquely identified by the color of the line. Moreover, 
colors provide an unambiguous way to glue the two simplices, using the unique gluing which respects all the induced 
colorings on the sub-simplices. This prescription associates a topology to any colored graph and it can be showed 
that all the graphs generated here are topological pseudo-manifolds in dimension d. 

Graphs of degree zero have the topology of the d-sphere. Moreover, their combinatorics can be described quite 
precisely [l8|. They are found as the graphs which maximize the number of faces at fixed number of vertices. For 
d > 2, there is a way to locally maximize the number of faces by inserting on any line a patch with two vertices 
connected by d lines. This contributes to d(d— l)/2 faces with exactly two vertices and does not change the degree of 
the graph. The family of closed melonic graphs is the family of graphs obtained by recursive insertions of this patch 
on any line starting from the closed graph with two vertices connected by (d + 1) lines. A typical melonic graph is 
depicted in Fig. [1] 

Equivalently, the generating function for melonic graphs satisfies some algebraic equation. Denote U the full 
connected 2-point function, (Ta^... ad f^\ bd ) = N~ d l 2 U8 13 Y\ c 8 aa b a i and similarly £ the 1-particle-irreducible 2-point 
function. One has as usual U = 1/(1 — £). In the melonic sector, one gets an additional equation to close the system 
on U, 




(2.13) 



This means that all contributions to £ factorize on d internal lines connecting the left and right vertices, with insertions 
of melonic 2-point graphs. One gets all contributions by inserting the full 2-point function itself U on each line. As a 
result, 

U = 1 +gU d+1 . (2.14) 

The solution to this equation is represented in the plane (g, U) in Fig. [5J Starting from the Gaussian model for 
g = 0, U = 1, one follows the curve of the solution until a stationary point is reached. There g — g c goes quadratically 
with (U — U c ) as the first derivative vanishes. Therefore (U — U c ) ~ ^Jg — g c . This singular behavior defines the 
thermodynamic limit and one extracts the entropy exponent 7 me ions = 1/2. 



C. Hard dimers on melonic triangulations via a tensor model 

The above model can be extended to include dimcr configurations dressing the dynamical lattices. A set of tensors 
(Xai—ad)t=o,...,d) together with their complex conjugate, is introduced to represent a dimer on a link of color (i). The 
hard-core repulsion is then fully encoded into the interaction of these tensors with the ones for pure lattices. As only 
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Figure 2. The equation U = 1 + gU d+1 in the plane (g, U) for d — 3. The Gaussian model is for g = 0, U = 1. The continuum 
limit is the black dot where dtf/dC/ = so that (U - U c ) ~ y^c - 3. One has f/ c = (d + l)/d and # c = d d /(d + l) d+1 . 



a single dimer can touch a site of the lattice, we need the following interaction 



£,dimcr(i) 



tax V y W Y\T U) 



(2.15) 



and similarly with the complex conjugate S_ 



r(t) 



In words, one removes the tensor of color (i) in the interaction 



Si and put the dimer tensor x^' °f the same color instead, and weights this vertex with y/gz. Dimcrs also need a 
quadratic part Sq™d' which is the same as Squad with \ instead of T. 
The partition function is 

Z = e- N " 3 ( z >ti = j [dT^d x {i) df^dx {i) ] exp -N d ' 2 (s quad + S^£f + S+ + S- + Sf nicr + 5 dimor ) , (2.16) 



where S± mer = Yli=o S± m ^ . We have denoted the free energy of the model 3(z, g) because it is indeed the grand- 
canonical partition function for hard dimers on the set of random bipartite (d+l)-colorcd graphs, with lattice chemical 
potential y/g and dimer activity z. This grand-canonical partition function has a degree expansion which is completely 
independent of the dimer configurations and hence exactly the same as for pure random lattices. 
In particular, at degree zero, 



So(z,fl) 



E 



g v ' 2 Z \°\. 



(2.17) 



mclonic graphs Gv 
dressed with D G T>(G\ 



III. CRITICAL BEHAVIORS 



Wc now somehow repeat the procedure of Sec. Ill Bt with details, to solve the model at the melonic order. Due to 
coloring, it is guaranteed that the full connected 2-point functions can be written 



/ T W fW \ _ N -d/2 rr 1 M 9 W \ _ N -d/2 w 

\J-ai—ad - L ai—a ( i) — Jv <~> , Uai-Od Aa r -a d / - " "j 

/ («) f(i) \ _ lrp(i) Ai) \ _ M-d/2 y 
\Xai---a d - L a 1 ---a d l — \- L a 1 ---a d \ai---a d / * ) 

all other components vanishing. 

First we use some Schwingcr-Dyson equations to relate the free energy to these 2-point functions. 



1 



Z N d 



[dT^d X (i) dT^dx (i) ] 



d 



ai,...,a d u ^ a 1 ---a d 



' \ ± ai - a d 



N d/2 (S quad (T,f )+S quad (X,X)+S+ +S- +Sf •»" +S di»") 



Since each term of the action is at most linear in each T®, evaluating the derivatives yields 

1 - U - N- d ' 2 (S+ + Sf mcr - ^ imor(l) ) = 0. 
We sum this equation with the same one coming from Tw and also sum over colors, 



N- d / 2 N~ d / 2 
1-U — (S+ + S- + S L l imci + 5 dimor ) + — — -(Sf mci + S d _ imer ) = 0. 



(3.18) 
(3.19) 



= 0. 
(3.20) 

(3.21) 
(3.22) 
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Then the connected expectation values in the above equation are recast as derivatives of the free energy to give 

(3.23) 



og a + 1 oz 



This reasoning is repeated from a Schwinger-Dyson equation with a x insertion instead to get 

z dE 



1-V- 



= 0. 



d+1 dz 

By combining the two above equations, we are able to write the derivatives of the free energy as 



dE 
dg 



z^f = (d + l)(l-V), 
oz 



which reduce the problem to finding U and V in the melonic sector. 



(3.24) 



(3.25) 



We denote £ 



full connected 2-point functions, 



A. The thermodynamic limit 

) the matrix of the 1PI 2-point functions. Its geometric series is exactly the matrix of the 



U W 
W V 



1 — Sj/ — £jy 

-Sit' i — £v 



1 
1 



Since they are symmetric matrices, one can just consider three real equations like 

U(l - Eu) - WT, W = 1, 
V(l - £y) - WY>w = 1, 
W(1-E V )-UY, W = 0. 



(3.26) 



(3.27) 
(3.28) 
(3.29) 



Therefore we just have to find some extra equations in the melonic sector which enable to close the system on 
U, V, W. The melonic sector is characterized by the fact that all contributions comes from 2-point insertions which 
are themselves melonic on any line. To generalize (|2.13l) in the presence of dimers, we just have to list all the possible 
ways of putting dimers in the graph of (|2.13l) , weighted appropriately with the dimer activity and some combinatorial 
factors counting the equivalent ways to put a dimer. This gives for T,jj 




dg\fz 




(3.30) 



gU d + 2dg^WU d - 1 + dgzVU d - 1 +d{d-l)gzW 2 U d ~ 2 , 



(3.31) 



where we have represented graphically the contributions at d = 3 with dimers as dashed lines. For Sy, one cannot 
make any dimer 2-point insertion so that there is only one kind of contributions, 



Sy = gz 




gzU d 



(3.32) 
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Figure 3. 



As and are independent of V, and is only linear in W, one can extract from (|3.29[) W in term of U, 

This provides an expression Ew(f) by inserting it into f|3 . 33|) . Then, one uses fl3.28[) to extract V in terms of U, 

1 — gz U a 

This way, one gets rid of V and W easily and reduces the algebraic system to a single equation, i.e. (|3.27[) where 
W, V are functions of U. It is useful to change U to u = zgU d . Then the equation becomes linear in g 1 ^, 

d+i i i du 3 + [(d + l) 2 z + \]u 2 - [2(d + i)z+ l]u + z 

z —3~ nd = U d _ __ . _ (3.36) 

y (du + l)((d+ l)u - l) 2 

The continuum limit is defined just like for pure random lattices as the points where u becomes singular, i.e. dg/du = 0. 
This is equivalent to the stationarity of g 1 ^ which reads 

dg\ „ (u-1) 2 +dhi 2 (u 2 -2u + 2) + d(u-l)(2u 2 -l) + d 2 u + d 2 u 2 (u 2 + u-3) , n 

— — =0 z = u o . (3.37) 

du ((d + l)u-l) (d(d-l)u-l) 

This gives the activity as a function of the parameter u. Plugging it back in ()3.36[) we also get the coupling g as a 
function of u. Therefore we obtain the line (z*(u), g*(u)) of the continuum limit as a parametrized curve. 



To get a picture of this line, one draws (|3.36p in the plane (g, U) for different values of the activity z, see Fig. 3(a) 



at d = 3 with z between —0.0935 and 4. Starting from z = one follows the local maxima obtained for different 
activities, represented as the red dashed line on Fig. 3(a) At fixed z, the first derivative with respect to U vanishes 
and the second does not. Hence, 

7 = 1/2, (3.38) 
for small enough activities. This is the behavior of pure random lattices, as expected when dimers are not critical. 



B. Multicritical behavior 



Let us focus here on d = 3. It can be observed on Fig. 3(a) that for large enough activities there is no local 
maximum anymore, so that the line of the continuum limit stops somewhere. To understand this, one should notice 
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that in addition to the set of local maxima which defines the physical continuum limit, there is a set of local, unphysical 
minima for each z, which actually come closer and closer to the physical line when z is varied. This can be seen on 
Fig. |3(b)| where it becomes clear that there exists a critical value z c of the activity where the physical line meets the 
unphysical set of local minima. At this point not only the first derivative of g with respect to U vanishes, but also 
the second derivative. Hence performing an expansion around this point at z = z c fixed we find 

g-g*(z c )~ A (U-U c ) 3 U sing ~ (g - , 9 *(^)) 1/3 7 = 2/3. (3.39) 

This means that the universality class of the continuum limit is not the one of pure lattices anymore. It is changed 
because dimers become critical at this activity z c . Beyond z c the solutions to dg/dU — move to the complex plane. 

Analytically this behavior take place when dz*/du = 0. Taking the derivative of (|3.37p and after removing an 
obvious root at u = —1/d we get 

1 + d + (-2 + d + 3d 2 ) u + (1 - 6d - 6d 2 + 3d 3 - 4d 4 ) u 2 + (5 + 2d + 2d 2 + 9d 3 - 2d 4 )du 3 

+ (-1 + d - 3d 2 - 5d 3 )du A + (d- l)(d + l) 2 d 2 u 5 = 0. (3.40) 

This is a polynomial of degree 5 for which we have not found any obvious root. Hence, they have to be evaluated 
numerically. At d = 3, the real roots are —0.075,0.185,1.37. The relevant one is the first, which corresponds to 
z r _ w —0.09, in agreement with the plot |3(b)| One also gets g*(z c ) sa 0.198, U*(z c ) w 1.61. Notice the critical act ivity 



is negative, as expected from the 2d case for non-unitary critical models. (For positive activities, one sees in Fig. 3(a 
that the solution of the model goes to g = 0.) 

Other interesting exponents are the ones which characterize the singularity when approaching it in the continuum 
limit. We move towards z c following the critical line (z*(w), g*(u)). Using some Taylor expansions it is found that z* 
and g* both start quadratically with (u — u c ), therefore 

g*{z) w g*(z c ) +b(z — z c ) 3 t 2 + regular terms. (3-41) 

To get the free energy in the continuum, we make use of the Schwinger-Dyson equations (|3.25|) and find 

dF dF dg„ dF , 

dz \z c oz dz og 

for some constant B rj —126.9. 



IV. MAPPING TO "COLOR-SENSITIVE HARD DIMERS" ON RANDOM BRANCHED POLYMERS 

It is known fill ] that planar graphs can be mapped to branched polymers, i.e. trees. Similarly melonic graphs are 
in one-to-one correspondence with a family of branched polymers. This is actually one way to understand the series 
of multicritical exponents of random matrix and tensor models which go like 7 = e — 1/m where e = 0, 1 and m > 2. 

While a tree is obtained from a planar graph by cutting lines, the mapping of melonic graphs is more abstract and 
goes as follows. First let us define a melon as a 1PI connected 2-point subgraph. At any scale, a portion of a melonic 
graph looks like Fig. i.e. a melon with left and right insertions on its external lines. One has two vertices connected 
by d lines, with all colors but, say, 0, and with (melonic) 2-point insertions on each line: this is a melon. One also 
has some (melonic) 2-point insertions on the "external" lines to the melon, with color 0. We map the two vertices to 
a single tree vertex and the lines of the graph to tree lines. The latter inherit the colors of the melonic lines. The 
tree vertex has a line of color which comes from its parent: it represents the line of color adjacent to the black 
vertex of the melon. Then the lines adjacent to the white vertex are mapped to lines connecting the tree vertex to its 
children. 

The process can be repeated. Indeed every melonic 2-point insertion appearing in Fig. U is itself of the type of 
Fig. [H This way one associates to any 2-point melonic graph a branched polymer with vertex degree (d + 2), and 
technically known as a connected, rooted (d + \)-ary tree. 

Now we can put dimers on melonic graphs and look at their mapping onto branched polymers. A dimer on the line 



of color 1 touching the white vertex on Fig. 5(a) is mapped to a dimer on the line of color 1 on the tree. As no other 
dimer can touch the white vertex, no other dimer can touch the corresponding tree vertex. It goes the same way for 
dimers on the lines of colors 2 and 3, so it might seem that hard dimers are mapped to hard dimers. However, this 
is not exactly true as we need to look at the effect of dimers on the lines of color 0. This case is shown on Fig. |5(b)| 
As one dimer touches the white vertex and another one touches the black vertex, they are mapped to dimers on the 
two tree lines of color which actually meet at the same tree vertex. 
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children 




Figure 4. The map from melonic graphs to branched polymers sends pairs of vertices to tree vertices. Every line of the graph 
is faithfully represented in the tree. The parent of a tree vertex comes from a melonic insertion connected to the black vertex, 
and children of the tree vertex represent melonic insertions connected to the white vertex. 



children children 




(a)A dimer, represented as a thick blue line, on the line of color (b)Two dimers on the external lines of a melon, here with color 
1 (or 2 or 3) is mapped to a dimer on the tree line of the same 0, are mapped to lines connected to the same tree vertex, 
color. In this case the hard-core repulsion is preserved. Therefore the dimers are hard-core except when they occupy 

lines of the same color. 



Figure 5. Dimers on melonic graphs are mapped to "color-sensitive hard-core dimers" on branched polymers. 



Therefore, the equivalent problem on branched polymer is a model of color-sensitive dimers. If a tree vertex has its 
parent of color i, then dimers on the children lines of color j ^ i have their hard-core repulsion preserved and cannot 
touch. But dimers carried by lines of the same color are allowed to touch each other. 

Then the 2-point function U calculated in the section IIIII is the generating function for these color-sensitive dimers 
on random, connected, rooted (d + l)-ary trees. Branched polymers become critical along the red dashed line of 



Fig. 3(a) which defines the thermodynamic limit in the grand-canonical ensemble of random branched polymers. A 



multicritical behavior is a point where trees and dimers become critical together, which happens at the critical activity 
z c as shown in Fig. |3(b)[ 



V. MORE DIMERS 



The parameter z has allowed us to observe some multicritical behavior. With more parameters, one can actually 
achieve higher order multicritical points as shown in [20^ . An interpretation of such behaviors has been given in p7| 
in terms of generalized dimers coupled to random lattices. However the free energy there was a generating function 
for melonic graphs coupled to dimers at fixed values of dV/2 — \D\, mixing the number of vertices and the number 
of dimers. This is not exactly what one would naturally consider: as the thermodynamic limit is defined by large 
values of the number of vertices V, one would like to work with a generating function which counts independently 
the number of vertices and the number of dimers. We show in this section how to write such a model, by a simple 
change of definition of the couplings in the model of f27j . 

The main idea of the model is to allow several dimers to touch if they connect the same two vertices, as shown 
in Fig. 6(a)| |6(b)| and to assign to a set of k dimers between two nodes some activity z^. These new terms arc 
introduced in the action as effective vertices which depend on the external lines (and their colors) which join the two 
nodes. For two nodes connected by d + 1 — k dimers, there are k external lines corresponding to tensors on the 
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(a)A graph dual to a triangulation at 

d = 3 (since nodes are four-valcnt) 
with a set of two dimers on the left, a 
single dimer in the center and a set of 
three dimers on the right. 



2 1 
(b)The dimer parts of the action can be drawn with internal 
(dashed) lines which represent dimers and internal faces (2 
lines and 1 face for the set of two dimers on the left, 3 lines 

and 3 faces for the set of three dimers on the right). The 
effective interaction only sees the external legs carrying some 
T and T. 



Figure 6. 



black vertex, for some choice of colors (cj), and k lines with the same colors for on the white vertex. This gives 



,c k } 



-gzd+i-kN' 



Jl-k)(d-k) 



tr Y[T<- C ' ] T<- 



Ci) 



(5.43) 



The index contractions are denoted with a trace, as it is otherwise a bit cumbersome. S 



{0,1} 
dimers 



aild inters are 



represented for d = 3 on the figure [6(b)] The terms with k = 1 can be seen as corrections to S qua d- The scaling with 
TV is chosen so that it is exact the same as gluing a node from S+ and a node from S- via d + 1 — k propagators, and 
hence (d + 1 — k)(d — k)/2 faces. The weight g comes from one y/g for the white vertex and one for the black vertex. 
The free energy of the model, at leading order, is a sum over melonic connected graphs 



^dimcrs(2,g) = ^2 g v ^2 II 

{G 2 v} D£V(G 2 v)k=2 



\D k \ 



(5.44) 



where \Dk\ is the number of sets of k dimers joining two nodes. In contrast with the model of [271 ]. g now really 
controls the number of vertices of the lattice. For simplicity, we assume that z± = 0, i.e. we do not allow the standard 
hard-core dimers considered in the previous section. 

Following the reasoning of [27|, similar to that of the above section, one has U = 1/(1 — E) and S can be obtained 
in the melonic sector as 



Y. = gU d 



d-2 
fe=0 



Zd-k U 



hence 



U 



d+l 



u - 1 



(5.45) 



(5.46) 



The continuum limit corresponds to the boundary of the analyticity domain of U in g where d(l/g)/dU = 0. This 
reads 



Zd 



dU d+1 -(d+l)U a + 



d(d-l)(d-2) 



d-2 



Z 2 U 



d - 1 + E[(^- 



-r 



k=l 



1 



Zd+i-k - (k + 1) 



Zd-k 



U K 



(5.47) 



This equation has the form <&(U, z^) = and implies that for generic values of the couplings it can be locally solved 
to get U in terms of the activities Z2, ■ ■ ■ , z&. Then the entropy exponent is the one of pure random melonic lattices 
7 = 1/2. 

However, (|5.470 is not invertible anymore if d^/dU = which is here equivalent to dzd/dU = 0. When this is true 
together with (|5.47p . the continuum limit becomes multicritical with 7 = 2/3 like in the previous section. Higher order 
multicritical points are also reachable when the first m derivatives of g vanish at the same time, d k (1 / g) / dU k = for 
k = 1, . . . ,m. Thus our conlusion is similar to that of [27| : this simple model inspired by dimers exhibits multicritical 
behaviors. 
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Figure 7. The plot of zz{U) which characterizes the continuum limit, for different values of Z2 = —10, . . . , 0. The black dot is 
the continuum limit at Z2 — z$ — 0. The red dashed line is the multicritical line where dz^/dU = and the black dotted line 
is the set of points where 1/g — 0. 



In [23], an example of a phase transition, i.e. a point where the physical solution to (|5.47[) meets another solution 
and that they exchange their roles, was presented for d = 6. Here we want to pay attention to d = 3 instead and 
compare with the hard dimers of the previous section. The equations for a multicritical point are 



1 _ U 4 + 3z 2 U 2 + z 3 U 
g~ U-l 



z 3 = 3 U 4 - 4 U 3 + 3z 2 U 2 - 6z 2 U, z 2 = -2U 2 , {5.4 



where we have ignored the possibility that U — 1. Clearly there can not be higher multicritical points as dz2/dll = 
would imply U = 0. Also notice that z 2 must be negative, in agreement with the expected non-unitarity of the dimer 
critical point. 

Fig. [7]plots z 3 as a function of U for different values of z 2 = —10, . . . , 0. The curves span a two-dimensional region 
of the plane. The black dot is the continuum limit at z 2 = z 3 = 0, where U = 4/3. The red dashed line represents the 
set of multicritical points where dz 3 /8U = 0, so in principle it can be reached smoothly from the initial configuration 
z 2 = z 3 = 0. 

However, there is another interesting curve which we have plotted as a black dotted curve and corresponds to points 
where the coupling 1/g vanishes, z 2 = — U 2 . It is not clear to us what the physical meaning of these points is, but it 
is clear that the multicritical line can not be reached from z 2 = z 3 = without crossing such a point where 1/g = 0. 



VI. CONCLUSION 



We have shown that hard dimers on dynamical triangulations, generated by random tensors in the large N limit, 
are solvable and exhibit a critical point in the thermodynamic limit. This is quite similar to the Yang-Lee singularity 
observed in the two-dimensional case, and in particular it also happens for a negative activity. The missing element 
of our analysis is obviously the identification of the universality class in term of a field theory. But we do not know 
yet the degrees of freedom in the continuum which emerge from the sum over melonic graphs. This is certainly an 
important challenge in the future. The present paper is a contribution to gain insights from the phenomenology of 
lattice models coupled to these melonic graphs. 

If one is able in the future to give a field theory description of the critical point, one could then consider a sort of KPZ 
correspondence to relate critical phenomena on flat space to critical phenomena coupled to geometric fluctuations. 



It is not clear how this correspondence could work in dimension higher than two. For example, it was shown in 25 1 
that the Ising model has no phase transition on dynamical melonic graphs while we do know it has one on flat space 
for any d. From this perspective a KPZ correspondence should be quite violent. However the present paper shows 
that non-unitary critical points can be observed, which makes the KPZ correspondence hypothesis more reasonable. 
Further work is needed to decide whether (non-)unitarity plays any role. 

Another challenge is to use tensor models beyond the melonic leading order. In matrix models, it is possible to do 
so using the Schwinger-Dyson equations for instance, to solve the higher genus contributions from the planar ones. 



For tensors, the algebra generated by the Schwinger-Dyson equations has been presented in [31| and contains an 
interesting sub-algebra generated by melonic observables, equivalently d-ary trees, that we propose to call the Gurau 
algebra 32j . It turns out that at leading order this sub-algebra reduces to (the positive part of) the Virasoro algebra 
as shown in [2(| . However this is not true at higher orders in the 1 /N expansion. The consequences of these different 
symmetries on critical behaviors should be studied too. 
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